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We study the quantum phase transition out of the Neel state in SU(3) and SU(4) generalizations 
of the Heisenberg anti-ferromagnet with sign problem free four spin coupling (so-called JQ model), 
by extensive quantum Monte Carlo simulations. We present evidence that the SU(3) and SU(4) 
order parameters and the SU(3) and SU(4) stiffness' go to zero continuously without any evidence 
for a first order transition. However, we also find considerable deviations from simple scaling laws 
for the stiffness and discuss scenarios to explain the observed behavior. 



I. INTRODUCTION 

The theory of quantum phase transitions^ has be- 
come a centerpiece in the research of quantum physics 
of strongly coupled condensed matter. It has found ap- 
plications in many branches of condensed matter physics, 
examples include: heavy fermion systems^, high temper- 
ature superconductivity^, quantum Hall effect, metal- 
insulator trans itio n, ultra-cold atomic gases^ and frus- 
trated magnetfP^. Aside from its many applications to 
natural systems, the study of quantum phase transitions 
has also turned into a full fledged theoretical endeavor 
in its own right, which although very mature in certain 
aspects^, is still in a stage of infancy when one considers 
the large number of different physical situations mod- 
ern condensed matter systems allow for. Indeed, most 
of the well understood quantum critical points (contin- 
uous quantum phase transitions) can be simply related 
to a classical universality class in one higher dimension. 
On the other hand, one of the most exciting directions 
in the field is the study of quantum critical points which 
lead to new universalities that are not natural to think of 
in the classical context. Such new criticality may arise, 
for instance, from essentially quantum phenomena such 
as the presence of low-energy fermions or complex Berry 
phases that have no natural classical analogue. 

An interesting theoretical proposal for a quantum crit- 
ical point that does not have a naive classical analogue 
was put forward a few years agcP^. The physical sys- 
tem is a two-dimensional quantum anti-ferromagnet on 
a square lattice with S=l/2 spins which transitions into 
a paramagnet with a broken translational symmetry, a 
valence bond solid (VBS). In a set of compelling argu- 
ments, it was shown that a naive application of the classi- 
cal theory which forbids a continuous transition could be 
invalidated by the presence of Berry phase terms, giving 
rise to a quantum critical point in a novel universality 
class. It is widely accepted, however, that a full non- 
perturbative understanding of the occurrence and prop- 
erties of the new quantum critical point are beyond ana- 
lytic reasoning and can only be established by unbiased 
numerical simulations. In a pioneering piece of worlP^, it 
was shown that these questions could be addressed in a 
sign problem free microscopic JQ model using quantum 
Monte Carlo techniques. Since then a number of works 
have studied this SU(2) symmetric modefSKS'. While all 



workers agree that the JQ model harbors a Neel state on 
one side and a VBS state on the other side, a clear pic- 
ture of the the nature and precise scaling at the quantum 
phase transition still does not exist. The most compre- 
hensive numerical study of the SU(2) JQ model however 
provides strong evidence in favor of a continuous tran- 
sition^, albeit with corrections to the naive scaling hy- 
pothesis. Corrections to naive scaling also appear in the 
study of the response of the JQ quantum critical point 
to single impurities^. Another studj^ for an unconven- 
tional quantum critical point in the spatially anisotropic 
bilinear bi-quadratic model has also found evidence for a 
continuous quantum critical point, though this model is 
less well studied. The only theory for this transition is 
also beyond a naive classical order parameter theorjff^. 
Excluding the JQ model and the anisotropic bilinear 
biqudratic 5 = 1 model, there is no known candidate 
for a quantum critical point in a numerically accesible 
two dimensional microscopic quantum Hamiltonian that 
does not map simply onto a higher dimensional classical 
field theory (the most well studied example of a quantum 
critical point with a higher dimensional classical mapping 
is the bilayer anti-ferromagnet, the mapping has been 
studied thoroughly^. Sorting out these issues is a very 
demanding one numerically, but is extremely important 
for the field of quantum criticality. 

In this paper we present the results of extensive QMC 
simulations on a larger TV extension of the original 
SU(N = 2) symmetric JQ model. The first study^ of 
these models reported T = QMC results using the va- 
lence bond basis method for the SU(3) and SU(4) models 
and found results consistent with a continuous transition 
with conventional scaling for quantities associated with 
the VBS and Neel order parameters. In this work we 
approach this model with the powerful finite-T stochas- 
tic series expansion (SSE) method 20 which allows access 
to finite-temperature, larger volumes and most impor- 
tantly for our purposes here, gives access to the stiffness 
and susceptibility. We find that these quantities have 
large corrections to scaling, just like those found in re- 
cent studies of the SU(2) symmetric JQ models. This 
result does not contradict the findings of conventional 
scaling in RefP^ since their results were focused on the 
order parameters. Indeed for the Binder ratio associated 
with the Neel order parameter we too find conventional 
scaling works very well (see Sec. III). The corrections we 
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find are restricted to the susceptibility and stiffness (de- 
fined below), which are susceptibilities of the generators 
of the SU(7V) symmetry and which cannot be accessed 
in valence bond basis method so far, which was used in 
the previous study. The major new result of this paper 
is a study of the possible scaling forms for these quanti- 
ties (i.e. the stiffness) in the SU(3) and SU(4) symmetric 
models. We also present a number of scenarios which 
could explain the observed numerical results. Larger-iV 
extensions of the SU(2) physics may have applications to 
cold atom systems, but are important in their own right 
as they allow one to approach the analytically well under- 
stood large— N limit 9 . Other numerical work has studied 
the phase transition between Neel-VBS by varying the N 
in SU(iV) continuously^!], without the introduction of a 
four spin Q term. 



II. MODEL 

We define our model on a square lattice, which can be 
conveniently separated into an A and B sub-lattice. We 
define N flavors of fermions, f a identically on each site 
(1 < a < N). An SU(iV) rotation will rotate among 
these N flavors with the following convention: 



/o -> fa - U*pfp 



(1) 



where the * implies complex conjugation. Now to pick a 
spin representation in terms of the f a fermions, we fix the 
number of fermions on each site, f a fa = 1 on the A sub- 
lattice and flffi = (N—l) on the B sub-lattice. This fixes 
the size of the local Hilbert space to be simply N on each 
site, on both the A and on the B sub-lattice. We define 
these N states with the following sign conventions and 
indicate their transformations under SU(7V) rotation, 



\a) A = /t|0) (\a) A - 
\a) B = f a \F) (|o) B 



)a) 
1)b) 



(2) 
(3) 



where |0) denotes the absence of any fermions and \F) 
is a fully filled site. The transformation properties im- 
ply that the N, A sub-lattice states transform in the 
fundamental representation of SU(iV) and the N, B sub- 
lattice states transform in the conjugate of the fundamen- 
tal representation. The transformations also imply that 
the state J2 a \ a )A\ct)B transforms as an SU(7V) singlet. 
This is the main reason we have chosen different repre- 
sentations on the A and B sub-lattices, i.e. to allow the 
formation of a two-site singlet, following previous work 
on anti-ferromagnets^HUl. 

To construct a spin model, we write down SU(iV) in- 
variant four fermion interactions that maintain the num- 
ber of fermions on each site; Explicitly, there are two 
such terms ft a fi a f}pfjp and r ia fipfj a f}p, the first one 
is just the identity operator in the projected space, while 
the second term which we call, is a projector onto an 
ij singlet (i is on A sub-lattice and j is on B sub-lattice). 



Note that the matrix elements of Pij are very simple, 

(ai/3i\Pij\a 2 p 2 ) = S^Btdazfa (4) 

Since they are always positive when non-zero, we can use 
these operators to construct sign problem free models. In 
particular the model defined by Hij — — ^- is the famil- 
iar Heisenberg model up to a constant, after identifying 
a = 1 with up(down) and a = 2 with down(up) on the 
A(B) sub-lattice. We define the SU(JV) "JQ" model with 
the following conventions. 



H. 



JQ 



N 



N 2 



(5) 



ijkl 



where the first sum is taken over all nearest neighbor 
bonds on the square lattice and the second term is taken 
over all elementary plaquettes of the square lattice with 
ij and kl being nearest neighbor bonds. Note that so de- 
fined, Hjq explicitly has no sign problem: all off-diagonal 
matrix elements are explicitly negative. 

We now turn to the observables of interest in our study 
here. There are N 2 — l traceless Hermitian matrices, X°j 
which generate the SU(./V) algebra. Of these, N — l can 
be chosen diagonal (so-called Cartan subalgebra). We 
choose them so that they satisfy: Tr[X a X b ] = S ab /2. We 
can then work out formulae for the "uniform magnetiza- 
tion" and "staggered magnetization" for the Cartan gen- 
erators in terms of the operator n a which measures which 
of the | a) states is occupied on each site {n a — f^f a on 
the A sub-lattice and n a = 1— fif a on the B sub-lattice). 



(G) 
(7) 



Note that the sign in the middle is opposite from what 
you would have expected for the familiar 5* = 1/2 case, 
because of the way our Hilbert space is defined on the 
A and B sub-lattice. The specific choices for the Cartan 
generators is detailed in Appendix [B] An outline of the 
algorithm and method is given in Appendix |A"} 



III. BINDER RATIO 

The first quantity which we study is the so-called 
Binder ratio. It is defined as, 
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(8) 



In addition to the Monte-Carlo averaging, the (...) imply 
averaging over the N — l Cartan generators. It is well 
known that the Binder ratio has a scaling dimension of 
zero and hence can be written in the following simple 
scaling form close to a quantum critical point, 



D 



(9) 



3 



5(7(3) 5(7(4) 




J/Q J/Q 



FIG. 1: A zoom-in of the crossing of the binder ratio, B, 
for SU(3) and SU(4) symmetric JQ models. Error bars were 
determined by bootstrapping the data. The solid lines are 
polynomial fits on the data and are plotted to guide the 
eye. The crossing point converges well for the larger system 
sizes. For L = 64, 96, 128, 192 the data crosses nicely within 
the estimated error bars, allowing very accurate brackets for 
the quantum critical points. (J/Q)l = (1.9905,1.9920) and 
(J/Q)* = (11.235, 11.255). 
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FIG. 2: Fluctuations of the spatial winding number (stiffness) 
for SU(3) and SU(4) symmetric JQ models. This is the same 
set of simulations as for the Binder ratio data shown in Fig. [l] 
(this is a wider view, compare a;— axes). Clearly there are sig- 
nificant deviations from the simple scaling behavior expected 
for a quantity of scaling dimension zero. 



diagonal generators. 



This means that if we fix the parameter LT (see Ap- 
pendix[C]), we assume everywhere in this work that z = 1, 
that B should become volume independent when g = 
in the scaling limit (i.e. for large L). The simplest way 
to find whether this is true is to plot B(g) for different 
values of L and look for a crossing of the various curves. 
Such a plot for both the SU(3) and SU(4) JQ model is 
shown in Fig. [I] We have collected data on system sizes 
ranging from L = 16 — 192. While the data is accurate 
enough to detect corrections to scaling for the smaller 
sizes the crossing converges very well as the system size 
is increased and the four biggest system sizes cross nicely 
within our estimated error bars. From this crossing we 
have an accurate estimate of the value of the critical cou- 
pling (see caption of Fig. [I]), which is consistent with the 
values quoted in Ref. I19l [thanks to the larger volumes 
simulated here, our brackets are about a factor of 2 times 
more accurate for both SU(3) and SU(4)]. 



IV. STIFFNESS 

Configurations of our SU(Af) model in the SSE method 
(see Appendix |A| can be thought of as a set of A" colors 
of non- intersecting closed loops. We can measure the 
winding number of each of the N colors of loops both in 
space and in time and associate a spatial and temporal 
current for each configuration, with each of the N — I 



■a _ X a w x 
J x aa * * a 

j a = X a W T 



(10) 

(ii) 



The fluctuations of these quantities averaged over MC 
sampling and the N — 1 diagonal generators (W%) = 

((jx) ) an( l (Wt) = (Ut) 2 )' are °f g rea t interest, since 
they too are expected to have no scaling dimension just 
like the Binder ratio, B, discussed earlier, because they 
are susceptibilities of the generators of a conservation 
law (S\J(N) symmetry). Hence, for the same reasons 
discussed in the case of B, they are expected to have 
a crossing point. Before turning to a numerical study of 
their crossing, we note that by linear response theory, the 
fluctuations in spatial winding number are simply related 
to the familiar spin stiffness (p s = / (3) and the fluc- 
tuation in the temporal winding number are related to 
the spin susceptibility (\u — (W^)/3/L 2 ), we shall refer 
to these interchangeably. 

Fig. [2] shows (W%)(g) for different L, for both the 
SU(3) and SU(4) model. While the data does cross, it 
does so very roughly. There are clearly large corrections 
to scaling behavior. Are these just corrections to scal- 
ing or is the leading scaling behavior itself affected? To 
make further progress in understanding the deviations 
from scaling, we study the crossing point of (W%) for the 
sizes L and 2L in Fig. [2] The crossing points have an 
x— intercept (the coupling J/Q) and a y— intercept (the 
value of (W//)). We study these intercepts as a func- 
tion of L with the hope of being able to extrapolate the 
large-L behaviors. In order to estimate an error bar for 
these crossing points, we carry out a somewhat labori- 
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ous bootstrap method, in which we generate synthetic 
data sets of the same length as our data (drawn by ran- 
domly sampling our own data) , carry out least squares fit 
to this data (with polynomials of order two), and then 
compute the crossing point in each simulated data set. 
The RMS of the extracted crossing point data is then a 
reasonable estimate of the error in determining the x— 
and y— intercepts of the crossing point. Fig. [3] shows the 
J/Q value for the crossing points as a function of 1/L. 
The data clearly converges for large L to a finite value. 
It is re-assuring that this value is completely consistent 
with the brackets for the critical point extracted from the 
Binder ratio data (shown as thick lines on the y— axis of 
Fig. [3]). Next, we turn to the y— intercept of the cross- 
ings. Here things are much harder to interpret because 
the data does not saturate even for the largest system 
sizes studied here. Of course, one could always argue 
that this is due to a finite size effect, but given the large 
volumes simulated here, and the absence of any sign of 
saturation in the y— intercept, we shall assume that (W%} 
diverges as the volume increases. We would like to note 
here that this by itself does not indicate a first order 
transition. At a first order transition (Wj) should grow 
linearly (because the stiffness itself remains finite instead 
of going to zero). We see in the insets of Fig. [4] that there 
is no evidence for a linear divergence. We find that for 
our data on the SU(3) and SU(4) models, a weak power 
law divergence seems to fit better than the logarithmic 
divergence found for SU(2)R see Figs, [i] and [5J 



DISCUSSION 
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FIG. 3: The a;— intercept (J/Q) of the crossing of L and 2L 
for the stiffness shown in Fig. [5J plotted versus 1/L. The 
errors bars are estimated by a bootstrap method detailed in 
the text. The error bars get smaller for large systems be- 
cause these curves are steeper; the a;— intercept of the cross- 
ing can hence be determined more precisely. At 1/L = 
the thick red line shows the bracket for the critical point 
from the analysis of the Binder ratio data (Fig. [TJ. Clearly 
the x— intercept of the crossing converges to a finite value, 
and this value is completely consistent with the Binder ra- 
tio crossing. The dashed line is a non-linear fit to the form 
g - A/L a , with g = 1.992, A = 3.20, a = 1.47 for SU(3) and 
g = 11.255, A = 51.3, a = 1.48. A motivation for this form is 
given in Appendix |D| 



It is interesting to understand the theoretical reasoning 
behind the numerical analysis presented here. While we 
cannot offer an unambiguous theoretical interpretation, 
we make some inferences on what could be taking place. 
First we summarize our conclusions from what we have 
observed, 

• There is essentially no difference in the qualita- 
tive behavior for the various quantities under study 
from SU(2) for both SU(3) and SU(4). Quanti- 
tative differences in quantities such as the criti- 
cal exponents must clearly exist as the symmetries 
change. 

• The Binder ratio crossing seems to follow the naive 
scaling expectation of a quantum critical point with 
one relevant direction and all other directions irrel- 
evant. There are corrections to scaling but they are 
of the conventional type, i.e., they do not appear 
to affect the leading scaling behavior and vanish for 
large system sizes. 

• The stiffness (and susceptibility) have large correc- 
tions to scaling. The crossing point converges on 
the a;— axis to the same critical coupling deduced 
from the Binder ratio, whereas the value of the Lp s 



itself (y— axis) has a sizable volume dependence. 
If we assume the observed behavior in Fig. [4] is 
the leading scaling behavior, we conclude that the 
winding number fluctuations diverges sub-linearly. 
The fact that the divergence is sub-linear is very 
crucial, since it implies that the stiffness (which is 
the winding number fluctuation divided by L) itself 
goes to zero as one approaches the critical point, al- 
beit not as 1/L like one would have expected from 
a naive scaling analysis. It is difficult numerically 
to distinguish between a log or a very weak power 
law divergence, as we possibly have here, since the 
y-axes in Fig. [4] do not change very significantly 
over the entire range of sizes studied. 

Perhaps, the simplest explanation for the stiffness mea- 
surement is that the volume dependance of the stiffness 
curve is a finite-size effect and that it would saturate as 
the system size is made even larger than that studied 
here. This would be consistent with a conventional scal- 
ing hypothesis for the critical points studied here. While 
we cannot rule this possibility out, we believe it is un- 
likely since the crossing in the Binder ratio seems to have 
converged well for the system sizes under study. If this 
explanation is to be made viable, one would need a the- 
oretical explanation for why the finite-size corrections to 
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FIG. 4: Testing a linear and logarithmic divergences in the 
y— intercept of the crossing points. We have made the "fits" 
using only the two biggest system sizes (insets fit to A + Bx 
and main panels A + B\n(x)). This is thus more a test of 
the model than a fit, hence we refer to it as a "fit". The 
inset shows a linear "fit" (straight line on a lin-lin plot) and 
the main graphs a log "fit" (straight line on a log-lin plot). 
The insets clearly show that a linear divergence is completely 
inconsistent with our data. The main graphs show that a log 
form fits a little better, but is clearly not as good as the power 
law (see Fig. [5| 



the stiffness are as large as they are. 

If we assume that the observed behavior is indeed the 
asymptotic behavior, we are faced with another vexing 
questions: What could result in this unusual behavior 
for the stiffness? This behavior must signal the break- 
down of one of the standard scaling assumptions made. 
One scenario which would result in the divergence of 
Lp s is if the scaling function depended non-analytically 
on an irrelevant operator. For instance, imagine that 
Lp s = figL 1 /", j^) with uj > and hence g u formally 
irrelevant. Normally one assumes the scaling function is 
analytic in its arguments and can hence set the second 
argument to zero. If it was not analytic and depended 
for instance like f(x, y) « ~ for small y, there would be 
a power law divergence in Lp s as observed here. One im- 
portant consequence of this scenario is that a power law 
divergence is more natural than a log divergence, since a 
log divergence (within this scenario) would require a sin- 
gular dependance of the scaling function and an exactly 
marginal operator in addition. Of course this scenario is 
completely speculative and if it is indeed true, we leave 
it for future work to understand what this irrelevant op- 
erator is and why it results in singular scaling functions 
for the stiffness. 

An interesting discussion is how the correction evolves 
from the SU(2), SU(3) to SU(4) critical points. Since the 
first theoretical work about such critical points^, it has 
been assumed that the properties of the critical point de- 



FIG. 5: Same as previous figure, but now a test of the pos- 
sibility of a power law divergence. Again the straight line 
drawn is made by "fitting" the form Ax B to the two largest L 
points. This model fits the SU(3) data extremely well down 
to the smallest system size! For the SU(4) data the fit is also 
excellent, if one attributes the deviation of the smallest sys- 
tem size to a finite-size effect. The unfamiliar numbers on the 
y — axis are because we have chosen to tick and label the log 
axis with points on a base of f .f , which although convenient 
here, is not one of the familiar bases. Doing a standard linear 
regression analysis on the entire data set for SU(3) and the all 
the data except for the L = 16 point for SU(4), we find esti- 
mates for the exponent, B 3 = 0.20(2) and B 4 = 0.16(2). The 
errors quoted here are based on roughly fitting the power law 
form to different data sets and making an estimate for how 
much the exponent differs. 



pend smoothly on N. Our findings are not inconsistent 
with that assumption. While it is numerically difficult to 
distinguish between a weak power law (as we might have 
here) and a logarithm, there are important theoretical 
distinctions that arise from such behavior. In particu- 
lar, if the log behavior found for SU(2jP31 results from a 
marginal operator, it is natural to assume that the scal- 
ing dimension of this operator should depend on N. This 
would make it hard to explain a log divergence in the 
SU(3) or SU(4) critical points, even though such behav- 
ior is not inconsistent with our numerical data. It could 
also be possible that in the SU(2) case a weak power law 
divergence of the same form as proposed here for SU(3) 
and SU(4) is present. A complete theoretical scenario is 
required to address these interesting possibilities. 
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FIG. 6: An identical analysis as that leading up to Fig. [5] 
but now for the fluctuations of the temporal winding number. 
The power law "fit" does not work quite as well here. Again 
the "fit" here is carried out by fitting the power law form 
Ax B to the two largest system sizes. The exponents we find 
here are B$ = 0.18(2) and B4 = 0.17(2), where the errors 
are determined roughly by fitting different sets of the data. 
Re-assuringly, they are in agreement with the values found 
for the stiffness. See also Fig. [7] for a comparison between the 
susceptibility and stiffness. 



Appendix B: Choice of Cartan algebra 

In this appendix we list our choices for the N — 1 diag- 
onal generators for N = 3, 4. There are different conven- 
tions with which these can be chosen. For SU(2) there is 
only one diagonal generator and the natural choice is: 



1/2 






-1/2 



for SU(3) we used: 




'1/(273) 
1/(273) 
-2/(273), 



Finally for SU(4) we used: 



/1/(2V2) 



V 





-1/(2V2) 

1/(272) 





-1/(272)/ 



Motrunich, and T. Senthil. All the simulation reported 
here were carried out on the BCX cluster at the Cen- 
ter for Computational Sciences at the University of Ken- 
tucky. 



Appendix A: Method and algorithm 

The method employed here is the stochastic series ex- 
pansion, which is well documented in the literature^. In 
the way we have defined the Hilbert space, one can think 
of the configuration space as a set of loops (the loops 
travel vertically up in time on the A sub-lattice and ver- 
tically down in time on the B sub-lattice) with a color 
assignment from one of N colors. For the updates, we 
have generalized the "deterministic" method of Ref. [501 
Sec. II D, to both the four site Q term and SU(N), this 
is possible since as noted in Sec. [TTJ because of the simple 
form of the matrix elements we are guaranteed that all 
configurations which appear at a given order of J and of 
Q, appear with the same weight. In this method there 
are two kinds of updates (1) [loop update] the N colors of 
a loop can be randomly assigned a new color and (2) [di- 
agonal update] swap of identity operators with diagonal 
operators. 



/-V(2n/2) 



V 

/-l/(2\/2) 



V 





1/(2^2) 

1/(272) 




-1/(2^2)/ 





-1/(272) 





1/(272) 
1/(272)/ 



All these choices are made so that they satisfy the 
standard normalization, Tr[X a X b ] = <5 a b/2. 



Appendix C: Ratio of size and temperature 

Assuming z — 1 scaling, a choice has to be made for 
the ratio of LT/c. Of course any fixed ratio will work for 
the scaling properties, but in order to best use the efforts 
of our simulations we would like to work as close the cubic 
limit as possible LT/c = 1. The problem is that we do 
not know what c is in terms of our couplings J and Q, a 
priori. To circumvent this problem, we pick the temper- 
ature so that the fluctuations of the winding number in 
space (the re-scaled spin stiffness) and the fluctuations 
in the winding number in time (the rescaled spin sus- 
ceptibility) are of the same order of magnitude. Indeed 
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FIG. 7: The ratio of the fluctuations of the temporal and 
spatial winding numbers evaluated at the crossing points of 
the curves for the susceptibility. From the saturation to a 
finite (non-zero, non-infinite) value of the ratio at large vol- 
ume shown here, it is clear that in the large volume limit the 
asymptotic behavior is identical for (W-?) and (W*). The ra- 
tio is not expected to go exactly to 1 away from the LT/c = 1 
point. The number is of order one, because of our choice of 
LT as discussed in Appendix |C| 



we would expect them to be identical at the isotropic 
point. This is achieved most simply by choosing a cou- 
pling close to the critical point and the choosing some 
appropriate system size (small is good enough, so that 
the correlation length is larger than the system size) and 
then adjusting the temperature till the fluctuations of the 
temporal and spatial winding number become approxi- 
mately equal. We found that working in units with Q = 1 
that the approximate equality of the temporal and spa- 
tial winding numbers was achieved for L//3 = 4 for SU(3) 
and for L//3 — 10 for the SU(4) case, we hence used these 
values for the simulations reported in this manuscript. 



Appendix D: Crossing point 

In this appendix we study how the crossing analysis 
is affected by a multiplicative correction to the standard 
scaling function. 

Imagine we had a scaling function that had the form, 
L a f(gL a ). Now lets ask for what value of g x it crosses. 
Assuming the function / is analytic in its argument we 



find close to the critical point, g a 



_ /(Q) 



for the 



crossing of curves of sizes L and AL. So, the value of the 
critical coupling does go to zero (the true critical point) 



for larger and larger L, but it already receives corrections 
without any corrections from irrelevant operators. The 
leading behavior of the quantity at the crossing point, 
described by such a scaling function is simply, L a f(0). 
In order to test the functional form of g x we have carried 
out a non-linear fit to the data on the stiffness as detailed 
in the caption of Fig. [3] We find v m 0.65(5) for both 
SU(3) and SU(4) which is consistent with the value of 
v reported in a previous study^l (the error quoted here 
for v has been determined very roughly by attempting a 
number of non-linear fits dropping different sets of data 
points). Interestingly the v in the previous worH^ was 
extracted completely from the data on the order parame- 
ters. This gives some credibility to the idea that there is a 
multiplicative correction to the naive scaling for the stiff- 
ness. The analysis is very similar for ln(L)f(gL a ), where 
again g x s» ^pr^Wo) ^ or iar S e an d the function would 
be described by, ln(L)/(0) at leading behavior. Our g x 
data would fit this model reasonably too, since the multi- 
plication of a log affects the numerical values very weakly. 



Appendix E: Susceptibility 

In this appendix we study the fluctuations of the tem- 
poral winding number (the re-scaled susceptibility). We 
carry out an identical analysis to that presented in the 
body of the text for (W x ). 

The susceptibility has an identical crossing behavior to 
that of the stiffness shown in Fig. [2] there is a drift in the 
crossing point in both the x— and y— intercepts. We have 
carried out an identical analysis of the crossing of L and 
2L and we find that the x— intercepts (shown in Fig.[3]for 
the stiffness), do not depend very much on whether we 
look at the crossing of the susceptibility or the stiffness. 
The y— intercept docs however, so we have made a sepa- 
rate plot for this quantity in Fig. [6j We find that unlike 
in the case of stiffness there are finite-size correction to 
the power law model. Yet the power of the fitted power 
law for the largest system sizes is in agreement with our 
findings for the susceptibility, as detailed in the inset. 
Finally, as an independent check that these two quan- 
tities diverge the same way as we approach the critical 
point, we have plotted the ratio of these two quantities at 
the crossing point value determined by the susceptibility 
crossing. This is shown in Fig. [7j As expected we find 
that the ratio of the quantities clearly goes to a constant 
at large volume, providing clear evidence that the two 
quantities diverge the same way, as one would expect for 
the emergent z = 1 scaling (we use the word emergent 
here, since the microscopic model is clearly not invariant 
under space-time rotations!). 
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